Measuring and Analysing Marginal Systemic Risk Contribution 
using CoVaR: A Copula Approach 



Brice Hakwa 

Fachbereich C - Mathematik - Stochastik, Bergische Universitdt Wuppertal 

Manfred Jager-Ambrozewicz 

Hochschule fur Technik und Wirtschaft Berlin - University of Applied Sciences HTW Berlin 

Barbara Riidiger 

Fachbereich C - Mathematik - Stochastik, Bergische Universitdt Wuppertal 



Abstract 

This paper is devoted to the quantification and analysis of the marginal risk contribution of 
a given single financial institution i to the risk of a financial system s. Our work expands 
on the CoVaR concept proposed by Adrian and Brunnermeier [5] as a tool for the mea- 
surement of marginal systemic risk contribution. We first give a mathematical definition of 
CoVaRa ■ Our definition improves the CoVaR concept by expressing CoVaRd as a 
function of a state / and of a given probability level a relative to i and s respectively. Based 

s\L^—l 

on copula theory we connect CoVaRa to the partial derivatives of Copula through their 
probabilistic interpretation (Conditional Probability). Using this we provide a closed for- 
mula for the calculation of CoVaRa L ~ l for a large class of (marginal) distributions and 
dependence structures (linear and non- linear). Our formula allows a better analysis of sys- 
temic risk using CoVaR in the sense that it allows us to define CoVaRa L ~ l depending 
on the marginal distributions of the losses L % and L s of i and s respectively on the one 
hand and the copula of L % and L s on the other hand. We discuss the implications of this 
in the context of the quantification and analysis of systemic risk contributions. We will, for 
example, highlight some of the effects of the marginal distribution F s of L s , the dependence 

parameter p, and the condition C (i l ) on CoVaRa^^ ^ ■ 
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1. Introduction 

With the last crisis it became clear that the failure of certain financial institutions (the 
so called system relevant financial institutions) can produce an adverse impact on whole 
financial system. The inability of standard risk-measurement tools like Value-at-Risk (VaR) 
to capture this systemic nature of risk (since their focus is on an institution in isolation: 
micro risk management) poses a new risk-management challenge to the financial regulators 
and academics. We can summarise this into two questions: 

1. How to identify System-relevant Financial Institutions ? 

2. How to quantify the marginal risk contribution of one single financial institute to the 
system ? 

As an academic response to this problems, Adrian and Brunnermeier proposed CoVaR ([2]) 
as a model to analyse the marginal adverse financial effect of a distressed single financial 
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institution on the financial system. They defined the risk measure CoVaR as the Value at 
Risk (VaR) of the financial system conditional to the state of the loss of a single institution 
and quantify the institution's marginal risk contribution (how much an institution adds to 
the risk of the system) by the measure ACoVaR. This is defined as the difference between 
CoVaR conditional to the institution under distress and the CoVaR conditional to the 
institution in a normal state. 

Thus the implementation of CoVaR involves variables characterising a single financial 
institution i (e.g. L l ) and the financial system s (e.g. L s ) respectively and variables char- 
acterising the interdependency structure within the financial system and between single 
financial institutions and the financial system s. This macro-dimension of CoVaR allows 
the integration of the dependence structure of i and s in the risk-measurement contrary to 
the standard risk measures ("micro- risk measure" e.g. VaR) where only variables character- 
ising the financial institution alone are considered. The CoVaR concept can be thus used 
by regulatory institutions as a macro-prudential tool (or as a basis for the development of 
other tools) to identify systemically relevant financial institutions and to set the adequate 
capital requirements. 

But its calculation represents an open problem. 

But its computation represents an open problem. Although some approaches have been 
proposed, Adrian and Brunnermeier [2] proposed for example an estimation method based 
on "linear quantile regression", Gauthier et al. [28] adopted a simulation based approach, 
Jager-Ambrozewicz [35] developed a closed formula for the special case that the joint dis- 
tribution of financial system characteristic variable is of the Gaussian type. In all these 
approaches there are some difficulties to flexibly model the stochastic behaviors of financial 
institution's specific variables and their dependence structure (interconnection) within a fi- 
nancial system, since only linear dependence are considered. 

Our aim is thus to provide a more flexible framework for the implementation of the CoVaR 
concept which allows the integration of stylised features of marginal losses as skewness, fat 
tails and interdependence properties like linear, non-linear and positive or negative tail de- 
pendence. To do this we first propose an improved definition of CoVaR which makes it 
mathematically tractable (see def. [3]) , and based on copula theory we propose a general an- 
alytical formula for CoVaR (see Theorem [4]). We use our formula to make some theoretical 
analyses and computations related to CoVaR. 

s \L* —I 

We conclude this article by applying our formula to compute the CoVaRa in the Gaus- 
sian copula (Section 5.1), t-copula (Section 5.2 1, and Gumbel copula setting (Section 5.3.1) 
respectively. We discuss the results of our computation and draw from this some interesting 

conclusions. We also give a general formula for CoVaRa L ~'in the Archimedian copula 
setting (Section 5.3). 



1.1. Definition of CoVaR and ACoVaR 



We recall here the definition of the value at risk (VaR) in order to define the CoVaR s][l 
as a conditional VaR following Adrian and Brunnermeier (0). 

Definition 1 (Value at Risk). Given some confidence level a € (0, 1) the VaR of a portfolio 
at the confidence level a is given by the smallest number I such that the probability that the 
loss L exceeds I is no larger than (1 — a). Formally 

VaR a := inf {I G M : Pr (L > I) < 1 - a} 
= inf {/ 6 R : Pr (L < I) > a}. 

In order to give a probabilistic interpretation of VaR a , we will employ the notation of 
quantiles as provided in the following definition (cf. [32] def. 2.12). 

Definition 2 (Generalised inverse and quantile function ). 

1. Given some increasing function T : R — > M., the generalised inverse of T is defined by 
T(y) := inf {a; ER:T(x)> y}. 

2. Given some distribution function F , the generalised inverse F^~ is called the quantile 
function of F. For a € (0, 1) we have 

q a (F) = F^ (a) := inf {x e R : F [x) > a} . 



2 



Note that, if F is continuous and strictly increasing, we simply have 

q a (F) = F- 1 (a) , (1) 

where F _1 is the (ordinary) inverse of F. Thus suppose that the distribution F of the loss 
L is continuous and strictly increasing. It follows 

VaR a = F- 1 {a) . (2) 

We note that typical values taken for a are 0.99 or 0.995. 

Assumption 1. Henceforth we consider only random variables which have strictly positive 
density function. Also in case we consider a bivariate joint distribution H (x,y) we assume 
that it has a density and its marginal distributions have strictly positive densities. 

So due to this assumption all considered distribution functions F are continuous and 
strictly increasing. Such an F is thus invertible and F~ x denotes the unique inverse of F. 
Let L % be the loss of the financial institution i and L s the loss of the system s without 
the institution i. At least since the financial crisis it is clear that the dependency between 
the system and the institution i must be analysed more seriously. A step towards such 
an analysis is done by explicitly defining CoVaR. Adrian and Brunnermeier denote by 

CoVaRa ^ the value of an institution s (or a financial system) conditional on some 

event C (L 4 ) depending on the loss L l of an institution i. Thus CoVaR^^ can be 
implicitly defined as the a — quantile of the conditional probability of the system's loss. 

Pr (l s < CoVaRT^^C (L l )\ = a. (3) 



They analysed in their work [5] the case that the condition C (L 4 ) refers to the loss L l of 
institution i being exactly at its value at risk or more generally being exactly at some specific 
value I. We have in this case in the context of ([3| the following expression, 



Pr (L s < CoVaRi v=l \V = l) = a. (4) 

Due to assumption [T] 

Pr {V = I) = 0, for any I G R. 

However we can define in the context of assumption [lj a conditional probability of the form: 
Pr (L s < h\V = I) for fixed I as a function of h as follows [cf. 14 p. 72) or[24 p. 71 ]. 



Pr (L s < h\L l = I) =: R t (h) 

, fi(l) 



dy. (5) 



Where /, (x) = f [x, y) dy is the marginal density of L l . 

Note that |5]) is defined only when (I) ^ 0; however, if S = {(l,y) '■ ft (I) ^ 0}, then 
Pr((L\L s )eS)=0. 

Remark 1. Due to assumption [T] we have that, 

• the functions Ri is well defined, (since fa (I) > 0, V I £ K), 

• Ri (h) is strictly increasing and continuous. 

As Ri (h) is strictly increasing, it follows that its is invertible. Based on this we provide a 

alternative definition for CoVaRa L ~ l which is more tractable from a mathematical point 
of view than that proposed by Adrian and Brunnermeier. 

Definition 3. Assume that L l and L s have density which satisfy assumption^ . Then for a 
given a G (0, 1) and for a fixed I, C'oVaRa L ~ l is defined as: 

CoVaR s J Lt = l := mf{h G M : Pr (L s > h\L l = I) < 1 - a} 
:= inf {h G M : Pr (L s < h\L l = l) > a} 

= R7 1 («) • (6) 
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Definition 4 (AC oV aR^ 1 ) . Adrian and Brunnermeier denote by ACoVaRa* the differ- 
ence between CoVaRa C ^ ^ condition on the institution i being under distress and the 
CoVaRa ^ condition on the institution having mean loss. 

ACoVaRf = CoVaR s ^ =VaR « - CoVaRf^^ . (7) 

ACoVaR sli is used as measure to quantify the marginal risk contribution of a single insti- 
tution i to the risk of the system. 

We will find in the next a closed analytical formula in terms of copula in the context of 
definition [3] (see Theorem [1]) . 

2. A Brief Introduction to Copulas 

In this section we introduce the notion of copula and give some basic definitions and 
important properties needed later. Our focus is on properties that will be helpful when 
connecting copulas to conditional probabilities and analyzing CoVaRa L ~ l and ACoVaR s \ l 
(for detailed analysis of copulas, we refer the reader to e.g. [36], [42], [44] or [47] and the 
references therein). 



2. 1 . Preliminary 

In order to introduce the concept of a copula, we recall some important remarks upon 
which it is built. 

Remark 2 (cf. |42j proposition. 5.2). 

1. Quantile transformation. If U ~ U (0, 1) is standard uniform distributed, then 

Pr (F -1 (Z7) <x)=F(x). 

2. Probability transformation. Assume F is a distribution function such that its 
inverse function F^ 1 is well defined. Let X be a random variable with distribution 
function F, then F(X) has a uniform standard distribution 

F(X) ~ U (0, 1) . 



2.2. Definition and basic properties of Copula 

Definition 5 (2-dimensional copula (cf. [H] def. 2.2.2)). A 2-dimensional copula is a 
(distribution) function C : [0, l] 2 — > [0, 1] with the following satisfying: 

• Boundary conditions: 

1) For every u € [0,1] : C(0,u) =C(u,0) = 0. 

2) For every u G [0, 1] : C (1, u) — u and C [u, 1) = u. 

• Monotonicity condition: 

3) For every (111,112) , (i>i,«2) € [0, 1] x [0, 1] with u\ < U2 and v\ < V2 we have 

C (u 2 ,v 2 ) - C {u 2 ,Vi) - C (ui, v 2 ) + C (ui, «i) > 0. 

Conditions (1) and (3) implies that the so defined 2-copula C is a bivariate joint distri- 
bution function (cf. [U] def. 2.3.2) and condition (2) implies that the copula C has standard 
uniform margins. We present now some important basic properties of copulas which we will 
use below (cf. 44J chap. 2). All this is summarised in the following theorem. 
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Theorem 1 (cf. [33] Thm. 2.2.7). Let C be a copula. For any v £ [0, 1], the partial derivative 
dC (u, v) I du exists for almost all u, and for such v and u 

Q< BCpv) <l 
du 

Similarly, for any u G [0, 1], the partial derivative dC (u, v) /dv exists for almost all v, and 
for such u and v 

Q< dCpv) <1 
dv 

Furthermore, the functions u i— >• dC (u,v) /dv and v H ► dC (u,v) /du are defined and non- 
decreasing everywhere on [0,1]. 

The following theorem makes the copula theory attractive as tool for stochastic modeling 
because it links joint distributions to one-dimensional marginal distributions. 

Theorem 2 (Sklar's theorem, cf. [Mj Thm. 2.3.3). Let H be a joint distribution function 
with marginal distribution functions F and G. Then there exists a copula C such that for 
all x, y E R U {-oo} U {+00} 

H(x,y) = C[F(x),G(y)}. (8) 

If F and G have density, then C is unique. Conversely, if C is a copula and F and G are 
distribution functions, then the function H defined by is a joint distribution function 
with margins F and G. 

This theorem is very important because it asserts that, using copula function, it is possi- 
ble to represent each bivariate distribution function as a function of univariate distribution 
function. Thus, we can use the copula to extract the dependence structure among the com- 
ponents X and Y of the vector (X,Y), independently of the marginal distribution F and 
G. This allows us to model the dependence structure and marginals separately. 

Remark 3. Assume (X, Y) is a bivariate random variables with copula C and joint dis- 
tribution H satisfying assumption [T] with marginals distribution function F and G. Then 
the transformed randoms variables U = F (X) and V = F (Y) have standard uniform 
distribution and C (U, V) is the joint distribution of (U, V). In fact 

C (u, v)=C {Pr (U<u), Pr (V < v)) . 

Corollary 3 (cf. e.g. [44] co. 2.3.7). Let H denote a bivariate distribution function with 
margins F and G satisfying assumption [7] . Then there exist a unique copula C such that 
for all (u,v) € [0, l] 2 it holds: 

C(u,v) = H(F- 1 (u),G- 1 (v)). 

3. Computing and Analysing systemic Risk Contribution with CoVaR^ L =l : A 
Copula Approach 

In this section we provide a copula based framework for the calculation and the theoretical 

s I — I 

analysis of CoVaRa as tool for the measurement of systemic risk contribution. To do 
this we will relate the notion of conditional probability to copulas and rewrite the implicit 
definition of CoVaR s \ L =l in terms of copula. Based on this we will derive some useful 
results. Specifically, we will obtain a closed formula which will provide a general framework 
for the flexible calculation and analysis of CoVaR s)[LZ=l in many stochastic settings. Based 
on this formula we will highlight some important properties of CoVaRa L ~ l and ACoVaR B ^ . 

3.1. Computation of CoVaRa L 1 using Copula 

s\L i —l 

We propose here in the following theorem a general framework for computing CoVaRa 
analytically. Our approach is based on the copula representation of conditional probability. 
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Theorem 4. Let L l and L s be two random variables representing the loss of the system i 
and institution s with marginal distribution functions F s and Fi respectively. Let H be the 
joint distribution of L z and L s with the corresponding bivariate copula C, i.e. 

H(x,y) = C(F i (x),F s (y)). 

Let us assume assumption^ and 

dC (u,v) 



g{v,u) 



du 



is invertible with respect to the parameter v. Then for all I € K CoVaR a 1 at level 
a, < a < 1 is given by 

CoVaR^^ 1 (a) = F' 1 (g- 1 (a, F t (/))) V a G [0, 1] . (9) 

s I — I 

Proof. Recall that the implicit definition of CoVaRa is given by: 
Pr (L s < CoVaR a ^= l \V = l) = a 
^Pr (F s (L s ) < F s (CoVaR a ^= l ) \F % (L l ) = Fi (I)) = a. 

Let V — F s {L s ) , U = Fi (L l ) , v = F s (CoVaR a lL ' =l ) and u = F, (I) i.e. 



Pr (L s < CoVaRt Ll=l \L l = V) = Pr (f s (L s ) < F s (CoVaR 8 ^ 1 ) \F t (V) = F, (I) 

= Pr(V < v\U = u) . 



Due to assumption[T]it follows from remark[2]that V and U are standard uniform distributed. 
In this case we can refer to ([TJ eq. (4.4)) and ([17] p. 263)) and compute the conditional 
probability Pr (V < v\U — u), as follows: 

Pr(V<v\U = u)= lim ^ (V < v,u < U < u + Au) 



A«^o+ Pr (u < U < u + Au) 

Pr (U < u + Au, V < v) - Pr (U < u, V < v) 

~ A™o+ Pr(U <u + Au) - Pr(U < u) 

C (u + Au, v) - C (u, v) 

= lim 

Au.-to+ Au 

dC{u,v) 

Now we are able to derive the explicit expressions of CoVaR a provided that, the function 
g is invertible with respect to the "non-conditioning" variable v. In this case we can write 
v as a function of a, u as follow 

v = g^ 1 ( a , u ) ■ 
Using v = F s (CoVaRi L ' =l \ and u = F t (I) we obtain 



F s (CoVaR 8 ^'- 1 ) = g- 1 (a, F< (I)) . 

Thus 

CoVaR s J L '= l = F- 1 (g- 1 (a, F { (I))) . U 

In practice the conditional level I for the financial institution i is implicitly defined by a 
given confidence level /3 such that 

I = F- 1 (fi) , (10) 



G 



P is specified by the regulatory institution. It represents the probability with which the 
financial institution i remains solvent over a given period of time horizon. Base on this 
information we can express CoVaR^- as follow: 

CoVaR^ L *= l = F~ x (g^ 1 (a, /?)) . (11) 
We remark that for a given marginal distribution of the system's losses F s the above ex- 

s\L^—l 

pression of CoVaR a has only as input parameter a and f3. This motivates the following 
definition. 

Definition 6. 

CoVaRi := CoVaR s a lL ^ 1 

s\L'=l 



Remark 4. Equation ( 11 ) is very important because it asserts that in the practice CoVaRo 
(or CoVaR^ ) contrary to standard risk-measurement tools like Value-at-Risk (VaR) does 
not depend of the marginal distribution Fi but depends only on the marginal distribution 
of the system's losses F s and the copula between the financial institution i and the financial 
system s. 

Remark 5. We can see from equation ^ that CoVaRa L ~ l is nothing other than a quantile 
of the loss distribution F s of the system s at the level a = g^ 1 (a, Ft (I)) i.e. 

CoVaR s } LZ=l = F- 1 (a) . (12) 



Equation ( 12 1 asserts that CoVaRa is just a value at risk of the whole financial system 
at a transformed level a = g^ 1 (a, F, ; (I)). This fact motivates the following corollary, which 

connects CoVaRa L ~ l to the value at risk at the level a of the financial system {VaR s a ). 
Recall that under assumption [l] the value at risk of the system at the level a of L s is in this 
case given by 

VaR% - F," 1 (5) . 

That is 

CoVaR^^ 1 = VaR s & 

Corollary 5. Provided that the function g (v, u) :— dC ^ v ^ is invertible with respect to the 
"non- conditioning" variable v, the equivalent confidence level a, which makes the Value at 

Risk of a financial system VaR s equivalent to the CoVaRu L ~ l at level a is given by: 

a = g^ 1 (a, u) with u — Fi (I) . (13) 

Hence, in general, given a condition quantile at the level a, we can find the corresponding 
unconditional quantile by transforming the conditional level a to a unconditional level a 

through the transformation function g^ 1 . Based on the fact that CoVaR s J L "'can be ex- 
pressed as a quantile. We can simplify the expression of CoVaRa L ~ l in a linear function 
when L s is assumed to have a univariate normal distribution. In fact if a random variable X 
follows a normal distribution with mean /x and standard deviation a. Then the transformed 
random variable Z = is standard normal distributed. This motivates the following 

proposition. 

Proposition 1. If the loss of the financial system L s is assumed to be normal distributed 
with mean fi s and standard deviation a s . Then 

CoVaR^'^ 1 = as®- 1 (a) + » s , (14) 



with a defined as in equation (13 1. Where $ denotes the standard normal distribution 
function. 

s\L^—l 

That means CoVaR a is in this case a linear function with respect to the transfor- 
mation <i> _1 (a). 



7 



Proof. Assume that L s is normal distributed with mean fi s and standard deviation a s . 
Let N s be the distribution function of L s then from ( 12 ) we have 



CoVaR s \ v=l = N' 1 (5) . 

And using the fact that any arbitrary normal distribution can be transformed to a standard 
normal distribution wc obtain 

CoVaR s \ L *= l = a^- 1 (a) + [m s ■ 

Corollary 6. Under the same conditions as the previous proposition, ACoVaRa 1 can be 
compute as follow 

ACoVaRf = a s (S" 1 (a*) - (a ro )) (15) 

Where and a m are the transformed level defining according to the corollary [5] when 
institution i is under distress and institution having mean loss respectively. 

In the following remark we summarise some properties of CoVaRa L ~ l as monetary measures 
of risk with particular attention to the concept of coherent risk measures [[6j [25] def. 4.5] . 
This summaries according to [B], properties that a good risk measure should have. 

Remark 6. As CoVaRA can be expressed as quantile of the distribution of the system's 

s I — I 

loss F s with respect to the transformed level a. It follows that CoVaRa as a function 
of a has the same properties like value at risk as a function of a level a. 

In particular following properties. 

Property 1. • CoVaRa L ~ l is a coherent measure of risk under elliptical distributions 
(cf. ex. 6.7). 

increases when the marginal distribution of the system (F s ) has leptokur- 
tosis (heavy-tailed) and positive skewness. (cf. ' L 5j IV. 2. 8.1). 

One important advantage of our formula is that, the expression of CoVaRa L ~ l (see 
eq. (|9|) can be separated into two distinct components. 

1. On the one hand the marginal distributions Fi and F s , which represent the purely 
univariate features of the single financial institution i and the financial system s re- 
spectively. 

2. On the other hand the function which represents the dependency structure be- 
tween the single financial institution i and the system s). 

This separation is very important for the analysis of systemic risk property of our formula. 
First, because it describes how the systemic contribution of one given financial institution de- 
pends on its marginal distributions Fi and the marginal of the financial system F s . Secondly, 
because it allows us to appreciate the effect of the copula of the systemic risk contribution. 

4. Tail Events and Systemic Crisis 

As asserted by Adrian and Brunnermeier [2] , the main idea of Systemic risk measurement 
is to capture the potential for the spreading of financial distress across institutions by gauging 
the increase in tail comovement (e.g. The prefix Co in [2] refers to conditional, contagion, 
or comovement). 

Definition 7. Forbes and Rigobon \27\] define contagion as a significant increase in cross- 
market linkages after a shock to one market (or a group of markets). 

Remark 7. During the crisis, the contagion effect appears to amplify the concentration of 
the financial system leading to an increase in probability that single financial institutions 
fail together with the whole financial system or that a large number of financial institutions 
fail simultaneously^ see Figure 1). 
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Figure 1: Bank Failures in the United States, from 1934 to 2009 (Source: |43] p. 443) 



The argumentation above highlights three important features of systemic risk. 

1. Systemic risk involves comovement. 

2. Systemic risk concerns a precise region of the involved losses distribution (e.g. the 
tail by Adrian and Brunnermeier [2: or the Distress region by 48J and 33J ) . Bernard 
et al. [9] and Hauptmann and Zagst [33] have taken this fact into consideration and 
proposed alternative definitions of CoVaR and systemic risk measure respectively. 

3. Systemic risk involves contagion. 

Hence in the context of the analysis and the measurement of systemic risk. The depen- 
dence between the financial institution i and the financial market s have to be considered 
only in a determined region of their joint distribution, (e.g. in the Tail or in Distress region). 

One way to do this would be to use dependence measures which allow the measurement 
of the dependence only in a defined region(cf. [40] chap. 6). For example in the case where 
the dependence structure is controlled by the correlation coefficient p, which is the case for 
elliptical copulas (e.g. Gaussian and t Copula). The conditional correlation coefficient has 
to be used instead of the unconditional correlation coefficient. 

Definition 8 (cf. |40| def. 6.2.1). Let U and V be two real random variables and A a 
subset of R such that Pr (V £ A) > 0. The conditional correlation coefficient p^ of U and 
V conditioned on V £ A is given by 



PA 



Cov (U,V\V e A) 



y/Var (U\V £ A) ■ Var (V\V £ A) 



Recall that the condition C in [2] refers to the loss being exactly at some specific 
value (e.g. VaRi) and because of assumption 1 we have that Pr (L l = ij = for any I £ K. 
To circumvent this problem we proceed as follows. Instead of considering the set where L l 
is assume some fixed value we follows Feller [21] p. 71 and consider the set where L l assumes 
values in an interval I = (1,1 + Al). We define 



P= 



lim 



Cov (L s ,L l \L l £ I) 



Ai^o y/Var (L s \L l £ I) ■ Var (U\L l £ I) 



(16) 



Remark 8. The use of the conditional correlation coefficient instead the unconditional al- 
lows the investigation of the effect of the Condition C on the systemic risk contribution. 

Another tool to measure the dependence of two random variables in one precise given 
region of their joint distribution is the so called quantile-quantile dependence measure A (a) 
introduced by cf. 17 . This is defined as 



A„ (a) =Pr(V> G- 1 (a) \U > F^ 1 (a)) 
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So according to the previous argumentation. The quantile-quantile dependence measure 
A (a) of L l and L s is per definition a natural indicator of the contagion between financial 
institutions over a threshold a (Note that, the typical value of a in our context are 0.99 or 
0.995 ). X u (a) = 0, for example, could mean that there is no contagion between i and s over 
the threshold a ). 

Remark 9. Let us consider a " lower- version" of the quantile-quantile dependence measure 
as defined by Roncalli [J7] (cf. 07] Remarque 58). 



A/ (a) = Pr(V < G- 1 (a) \U < F' 1 (a)) 



If we redefine the condition C (L l ) in the implicit definition of CoVaR,l C * (see. eq.([3]) 
and Q) by replacing "=" by "<" we have the interesting relation. 

Pr (l s < CoVaRf^lC (V)) = Pr ( L s < CoVaRf^ \L l < I 



"Assume that CoVaR^^ = F' 1 (a)" = Pr (L s < F' 1 (a) \V < Fr 1 (/3)) 
"Assume that a = P then we have " = Pr (L s < F^ 1 (J3) \V < Fr 1 (/?)) 

= A, (p) . 

The asymptotic consideration of X u (a) and A; (a) leads to the following definitions. 

Definition 9 (cf. f42| def. 5.30). Let [U, V) be a bivariate random variable with marginal 
distribution functions F and G, respectively. The upper tail dependence coefficient of U and 
V is the limit (if it exists) of the conditional probability that V is greater than the 100a — th 
percentile of G given that U is greater than the 100a — th percentile of F as a approaches 
1, i.e. 

X u := lim X u (a)= lim Pr (V > G _1 (a) \U > F^ 1 (a)) 

a— ¥l~ a— >1 

If A« G (0, 1] then (U, V) is said to show upper tail dependence or extremal dependence in 
the upper tail; if X u = 0, they are asymptotically independent in the upper tail. 
Similarly, the lower tail dependence coefficient A; is the limit ( if it exists ) of the conditional 
probability that V is less than or equal to the 100a — th percentile of G given that U is less 
than or equal to the 100a — th percentile of F as a approaches 0, i.e. 

X r .= lim A/(a)= lim Pr (V < G~ x (a) \U < F^ 1 (a)) . 

a— >0+ a— >0+ 

If [L 1 , L s ) does not show tail dependence (upper and lower) the extreme events of L l and 
L s appear to occur independently in each margin. This means that they are no-contagion 
betwenn i and s. 

Let us consider the bivariate Gaussian copula as model for (L J ,L S ). One can show that 
the bivariate Gaussian copula does not have upper tail dependence when the corresponding 



correlation coefficient p is smaller than one (see 17). As can be seen in Figure 2, regardless 
of how high a correlation we choose, if we go far enough into the tail, extreme events appear 
to occur independently in U and L s . That means the Gaussian copula is related to the 
independence in the tail. Thus Gaussian copula is not a good model for the analysis of 
systemic risk contribution between i and s. This is the reason why we connect the CoVaR 
concept to copula in order to develop a closed formula for CoVaR s][L =l allowing the analysis 
and the computation of systemic risk contribution for a more general stochastic setting than 
only the bivariate Gaussian setting as already done in (|35|). 



5. Applications 

In this section we apply the result developed in section [3] to compute and analyse 
CoVaRT and ACoVaR in some probabilistic settings. We first consider a general case 
where the joint behavior of L 1 and L s is modeled by a bivariate Gaussian copula. In partic- 
ular we will analyse here the case where the margins L l and L s are assumed to be univariate 
normal distributed. This special case (Gauss copula and Gaussian margnis) was already 
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considered in 35 but in a different approach. Jager-Ambrozewicz [35] assumes that the 
random vector (L l L s ) follows a bivariate Gaussian distribution. Then, based on the prop- 
erties of the conditional bivariate Gaussian distribution (cf. e.g. [Mj eq. 2.6), he develops a 

closed formula for CoVaRa L ~ VaR \ This approach imposes thus the univariate normality 
of both margins (L l and L s ). The method provided in this article is more flexible because 
it allows each margins independently of other to take a large class of distributions functions 
(for example we ca assume that L l is normal distributed and that L s is i-distributed). One 
other restriction of the formula proposed in [35^ and also the method presented in [2] is 
that, both do not take into account tail events and tail comovements since the Gaussian is 
asymptotically independent in both tails i.e. X u = A/ = . In fact we have (cf. [55] p. 17). 



X u = 2 lim 



2 lim 



l - $ 



l - $ 



a — pa 
'ay/1 - p N 



(17) 



from which it follows that 



if p < 1 

1 if p = 1. 



This presents a big gap since both phenomenons (tail events and tail comovements) are 




Figure 2: X u (a) for bivariate Gaussian Copula 

supposed to be the main features of systemic crisis (cf. [51]). Our formula covers this gap by 
allowing us to consider other dependence models, especially those which are appropriate for 
the modeling of the simultaneous (tail) behavior of losses during a financial crisis. So we will 
also consider the case where the dependence between L l and L s is modeled by a i-Copula, 
Gumbel-Copula. At the end of this section we will describe how to develop a closed formula 
for the computation of CoVaRa L ~ l for Archimedean copula. 

5.1. Computation of CoVaRa L ~ l in a Gaussian Copula Setting 

We assume here that the interdependence structure between U and L s is describe by 
a bivariate Gaussian copula. The bivariate Gaussian copula is defined as follows (cf. [H] 
eq. 2.3.6 ): 

C p {u,v) = ^ (®{u)- 1 ,§{v)- 1 ) , 

where $2 denotes the bivariate standard normal distribution with linear correlation coeffi- 
cient p, and $ the univariate standard normal distribution. 
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Hence, 

f*- 1 ^ r^W 1 / 2pst- 3 2 -t 2 \ 
G = / / ; exp — — dsat. 

Note that C p (u, v) can be express as 

In fact let X = (U, V) a standard Gaussian random vector with correlation p. Then we 
have: 

$2 (u, u)=Pr(f/<M,K< u) 

^ 1 f2pst-s 2 -t 2 \ , , 
-exp \ — ^ — ) dsdt 

-oo «/ — oo ' 1 



27r/I^p2- -V 2(l-p 2 ) 



this implies that 



a$ 2 (u, v) r 1 / 2put -u 2 -s 2 \ 

1 / — (s — up) 2 + p 2 u 2 — u 2 \ 

-oo WW^" ( 2(W*) J *» 

1 ^ -( S -u P f-u 2 (l-p 2 ) \ 

1 ( -u 2 -{s-upf\ , 
= exp | — — h ^— as 



exp ( ) / — exp ( — -q ] c?s 



t> — up 



= (/>(«)■$ 



where denotes the density of the standard univariate normal distribution. Therefore we 
have, 

$ a («,«) = J <t>(x) •* f-^==J d:E ' 
The expression of the bivariate Gaussian copula is then 

By making the substitution t = $ (x), we obtain 



By considering the expression |18|) we have: 
<9C P (u, v) 



9 («, w) 



g ( r" 1( ^$' 1(t) ) * 



On 



-.f*-^)-/*- 1 ^ (19) 
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contour of v=rf ' (u,q) with p = □ 
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0.2 0.4 0.6 0.B 




Figure 3: g 1 of the Bivariate Gaussian Copula 



The function g (v,u) is strictly monotone with respect to v. To compute its inverse, we set 
g (v 7 u) = a and solve for v. 

v = g- 1 (a, u) = $ (p<S>~ 1 (u) + ^/l - p 2 <S>~ 1 (a)) . 



Remark 10. If we set in the above equation p — 0, we obtain v = a for all u £ [0, 1] (see 



Figure 5.1 )). This is not a surprise because zero correlation means independence under the 
normal copula setting. 

So according to theorem [4] and the development make in section 4 we have the following 

s\ —l 

formula for CoVaRd when the dependence is modeled by a Gaussian copula. 



Proposition 2. Assume that the copula of L l and L s is the Gaussian copula then 



CoVaR s \ v = l = F' 1 ($ (p = $- x (Fi (I)) + ^l-pl®' 1 (a) 
In the context of remark^ we have 

a = $ fp=$~ 1 (Fi (I)) + ^Jl-pl®- 1 (a) 



(20) 



Where F l and F s represent the univariate distribution function of L l and L s respectively. 

Let us suppose as a particular case that L z and L s are Gaussian, that is F s = N s and 
Fi — Ni , where and N s are Gaussian distributions of the losses L % and L s with expected 
values pi, p s and standard deviation a s , Oi (This correspond to the case considered in |35|). 
This h We obtain the following closed analytical expression of C oV aR s)[L =l hi the Gaussian 
setting (Gaussian Copula and Gaussian Margins) 



CoVaR s \ L *= l = p= ^(l- Mi ) + ^l-pla^- 1 (a) + p s . 
In fact we have 

CoVaR s \ v=l = N- 1 ($ (p^' 1 {Ni (I)) + y/l - pl^ 1 (a))) 

= N- 1 (n s (osp^- 1 (N (/)) + a^l-pl^ 1 (a) + 
= o-sp^- 1 (Ni (I)) + a^l-pl^ 1 (a) + p s 

'"'"^ •-a s ^l^pi<S>- 1 (a)+p s 



(21) 



cr s p = $ 1 | $ 



I - Pi 
o-i 



",!>■ \ ) + (W 1 -fit® 1 ( a ) + Ms 

= P=— (I - Pi) + V 1 - Plos^ 1 (a) + p s - 
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And we have in this case (cf. [3^ eq. 6.1) 



p= = lim p — (22) 

A(^0+ / 2 . q V°r(L') 

yP -T (,1 P ) V ar(L*\Li£l) 

Remark 11. Note that p— can be either greater or smaller than p since Var (U\L l G /) can 
be either greater or smaller than Var (L l ). This fact have to be considered when analysing 
the effect of the dependence parameter p on CoVaR. The consideration of p— highlights 
also the impact that the change in volatility can have on the systemic risk contribution. 
This a very important since the behavior of the volatility are not the same depending if we 
are in distress region or not. In fact one can observe that distress times are characterised 
by high volatility. 

Corollary 7. Assume that L s and L l are Gaussian distributed and centered at zero then. 
CoVaR s ^= l = ( p = — J I + y/l - plcr^- 1 {a) 



= ( p=— J VaRi + y/l - plVaR s a whenl = VaR 



■p 



Remark 12. If p= = then CoVaR 3 ^ L ' =VaR " = VaR%. 

Let / be the value at risk of the single institution at the level j3 i.e. I — VaR l g = Pj (f3). 
Then we have the following expression of CoVaRi ( see def. [6|. 

Corollary 8. In the Gaussian setting, we have 

CoVaRi = p =( T s $- 1 (/3) + y/l-pLa,®- 1 (a) + p s . (23) 

Proof. 



CoVaRi = h ■ VaR l g - h ■ p t + y/l - picr^" 1 (a) + p s 



P= 



P= 



a, 



-as®- 1 ^) + ^1-pla^- 1 (a) + p s 



Remark 13. We remark that unlike in equation 21 the expression of CoVaR^, does not 
depend of the loss distribution's characteristic (e.g. standard deviation o~i and mean pi) of 
the financial institution i. 

Corollary 9. In the Gaussian setting. The map 

(a, /3) i — > CoVaR^ 
is increasing with respect to its both parameters. 

Now we refer to definition [4] to compute ACoVaRa 1 ■ The result of our computation is 
provided in the following proposition. 

Proposition 3. In the Gaussian setting, ACoVaRa* is given by 

ACoVaRf = p = a s <S> (py 1 . (24) 
Proof. According to definition [4] we have 

ACoVaRi = CoVaR a lL * =VaR > 3 - CoVaRf}^^* 



= p=— • VaR l B - P= — ■ pi + y/l - p 2 = VaR s a 
= p= — • VaRg - p = — ■ pi 
= p= — ■ (VaRl - m) 



0~s 0~ s 
P=— ' Pi ~ P= — 



P 

o-i p 

p=. — ■ (cr^" 1 {P)+fH - pi) 
o-i 

p^a^(p)- 1 . 
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Remark 14. From equation (24) we observe that if the financial institution i and the 



financial system s are not correlated, the risk contribution of i to s is zero. 

Let us impose now that, the loss of the financial system L s alone follows normal univariate 
distribution. Then according to proposition [T] we have 

CoVaR s ^ L '= l = a^- 1 (&) + fx s . 

Additionally if we also assume that L l is normal distributed. Then we have 

a = $ 1 (Ni (/)) + y/l -pi®- 1 (a)^ 



In sum we have 

CoVaR^- 1 = a s (V 1 ($ (p^ 1 (N, (I)) + ^T^pl®- 1 ( 
= a s (p=$- 1 (Ni (/)) + Vl-pi®- 1 (a)) + p s 

= <r s (p^- 1 U {^f) ) + V^Pl^- 1 («)) + A* 
2 - Mi 



a s p = P— ^ + v/Wi*- 1 (a) + 



Ps 

= p=—(l- Mi) + ^Vl-P^" 1 (") + Ms 

and 

CoVaRi = a s (V 1 ($ (/9=$^ 1 (/3) + v 7 ! - pl^ 1 (a) 
= a s (p=$- 1 08) + Vl-Pl^ 1 (")) + Ms 
= a,^*- 1 (0) + d^l-pi^ 1 (a) + Ms- 

Similary we can compute ACoVaRa 1 as follows. Recall (see corollary [sj) 
ACoVaRf = <j s (dd) - (cf m )) . 

And we have 



a m = $ [p^- 1 (0.5) + Vl-Pl^ 1 (c 

= $ f ^/wlr 1 (a)^ 



and 



a d = $ (p = $- 1 08) + v^r 1 (a) 

Hence 

ACoVaRf = (7 s (V 1 (* (p=$- 1 (/3) + v/l -pi®- 1 (a))) - (* (^/l - p| $-i (a) 

= a s ((p^- 1 08) + v/T^f" 1 (a)) - (v'Wlr 1 (a))) 
- ^p^" 1 09) • 

5.2. t-copula 

The Student t copula represents a generalization of the normal copula by allowing for 
tail-dependence through the degrees of freedom parameter. 
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Definition 10 (bivariate t distribution). The distribution function of a bivariate t- distributed 
random variable with correlation coefficient p is given by: 



tp. M (u,v) 



oo J — oo 



s 2 + t 2 - 2p t st 
^(1-(P) 2 ) 



dsdt, 



where v is the number of degrees of freedom. 

The Student t copula can be consider as a generalization of the Gaussian copula. He 
has in addition to correlation coefficient p a second dependence parameter, the degree of 
freedom v, controls the heaviness of the tails. For v < 3, the variance does not exist and 
for v < 5, the fourth moment does not exist. The t copula and the the Gaussian copula are 
close to each other in their central part, and become closer and closer in their tail only when 
v increases. Especially for both copulas are almost identic when v — > oo freedom parameter. 

Definition 11. The bivariate t copula, C l p v is defined as 
Cp,„ («,«)=*„,„ (t- 1 («) ,t" x («)) 

r^V) ft^M | / s 2 + t 2 



2p t st 



^(l-(p) 2 ) 



The tail-dependence coefficients the t Copula is given by (cf. e.g. 
of the symmetric property of t distribution we have, 



dsdt. 



eq. (5.31)) Because 



A; — A u — 2 — 2tj,_|_i 
From which it follows that, 



1 + P 



An — 



>0 i£p>-l 



if p 



-1 



Provided that p > 1. The bivariate t copula is thus able to capture the dependence of 
extreme values. 




□ 6 OB 




Figure 4: Tail Dependence Coefficient for t-Student Copula 

Following Roncalli [47] (cf. e.g. j47j p. 299) , we can express the t copula C* „ (u, v) as 
follows: 



Jo 



v + 1 



1/2 



v+ [tu 1 (u)]' 



t-Hv)-pt^(t ) 



dt. 



(25) 
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Now based on theorem [4] we compute the expression of g (v,u). We obtain 
g{v,u) 



dC l,u ( U ' V ) 



du 



u+1 



1/2 



On 

v+l \ V2 K 1 (v) - pt- 1 (u) 

,2 I / ~ 



The function g is invertible and its inverse is obtained by solving the equation 

\ 1/2 



g(v,u) = t v+x 



v + 1 



(«) - pK 1 («) 



^+ [i., 1 (u)] / \]l- {pf 



for v. This leads to, 



u = g 1 (a,u) = i„ 



pt v 1 («) + \ 



l-(p) 2 )(,+ [^( M )] 2 ) 

^ W («) 



V 



We obtain the following formula for CoVaRa L 1 and CoV^ai?^ when the dependence is 
modeling by a t-copula. 



Proposition 4. 



and 



( 



P^- 1 (Ft (I)) + 



V 


\ 






t 


( 


CoVaRi = F- 1 








V 


V 



\ 



1 - 



{p=f)(u+[t^{F l {l))} 2 ) 



\ 



v + 1 



J 



v + 1 K+1 [ ' 



J 



Where F l and F s represent the univariate distribution function of U and L s respectively. 
(3 denotes the regulatory risk level of the financial institution i. 

If L l and L s are t distributed with degrees of freedom v then 



CoVaRi L ^ 1 = (p= • 1) + \ 
and we have in this case (cf. [40] eq.6.B.31) 



(l-{p=) 2 ){v + P) 1 



(26) 



p- = lim f 5 r , I = (l:l + Al). 



P 2 



Note that the dependence in the Gaussian and t-copulas setting are essentially deter- 
mined by the correlation coefficient p (elliptical copula). The correlation coefficient is often 
considered as being a poor tool for describing dependence when the margins are non-normal 
(cf. [H] . This motivates the next section. 
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5.3. Bivariate Archimedean Copulas 

s\L* —I 

We can give in this section a general expression of the CoVaRT for some Archimedean 
Copulas. Archimedean copulas are often used in practice because of their analytical prop- 
erty, and ability to reproduce a large spectrum of dependence structures. Differently from 
the elliptical copulas, The definition of a bivariate copula are not derived from a given bi- 
variate distribution. The construction of Archimedean copulas is based on special function 
(the so called generator). The generator of a Archimedean copula is a convex and strictly 
decreasing continuous function tp from [01] to [0,oo] with tp (1) = 0. 

Definition 12 (pseudo-inverse, cf. [42] def. 5.41). define a pseudo-inverse of tp with domain 
[0, oo] by 



(t) = 



tp~ l {t) ifO<t<ip(Q) 
iftp(0)<t<oo. 



Note that the composition of the pseudo-inverse with the generator gives the identity i.e. 

ip^{if(t)) = t. vte[o,oo]. 

If <p (0) = oo the generator is said to be strict and it is equivalent to the ordinary 
functional inverse tp^ 1 . 

Given a generator tp we can construct the corresponding Archimedean copula as follows 

C(u, v) = tp [ ~ 1] (<p(u) + tp(v)) . 

The lower and upper tail dependence coefficient of an Archimedean copula can be com- 
puted using following corollary. 

Corollary 10 ( [44] co. 5.4.3). Let Let C be an Archimedean copula with a continuous, 
strictly, decreasing and convex generator tp. Then 

l-<p~ 1 (2x) 

A„ = 2 - lim 



A; = lim 



x-yQ+ 1 — tp 1 (x) 

1 - tp- 1 (2x) 

oo 1 — tp- 1 (x) 



Proposition 5. Let C be an Archimedean copula with a continuous, strictly, decreasing and 
convex generator tp i.e. 

C(u, v) = p- 1 (tp(u) + p(v)) . 

Then the function g defined as in theorem^ is given by (cf. \44j Thm. 4-3.8): 

= dC{u,v) = tp' (u) 

9[VlU > du tp>{tp-i[tp{u)+tp{v)}y 

Set g (v, u) = a and solver for v, we obtain the inverse of g. Namely: 



-x, x -iff ,_i/V(«m 

(a,u) = tp \tp\tp \—^))^ L P 



(«) 



Based on theorem [4] we derive the following proposition , which gives the expression of 
CoVaR 3 ^'^ for some Archimedean copulas 

Proposition 6. Let Let C be an Archimedean copula with a continuous, strictly, decreasing 
and convex generator tp Let L s and L l be two random variables representing the loss of 
the system s and institution i with joint distribution defined by a bivariate copula C with 
marginal distribution functions F s and Fi respectively i.e. 

F L i L s (x,y)=C(F i (x),F s (y)). 
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If C is an Archimedean copula with a continuous, strictly, decreasing and convex generator 



ip, then the explicit (or closed) formula for the CoVaRd at level a, < a < 1 for a 
certain fixed value I of L l is given by: 



CoVaR^ 1 = F- 1 (g- 1 (a,F t (l))) ■ 



ip \<p\<p 



,-x f<P'(Fi(l)) 



and 



CoVaRt = F7 l 



We are particular interested here by the archimedean copulas showing positive upper tail 
dependence(e.g. Gumbel copula). 

5.3.1. Gumbel copula 

Definition 13. The bivariate Gumbel Copula function is given by (cf. ex. 4- 25) 



Cf u (u,v) — exp I — (—lnu) n + (—lnv) c 



1< 9 < oo, 



where 9 represents the strength of dependence. Note that: 

• For 9 — 1 we have no dependency copula, i.e. Cf u (u, v) = uv 

• For 9 — > oo we have the perfect dependence i.e. Cg (u, v) = min (u, v) = M (u, v) 
with m and M represented the Frchet-Hoeffding lower and upper bound respectively. 

The generator of the bivariate Gumbel is given by ipg (t) = [—Int) for 9 > 1. 
The tail dependence coefficient of the Gumbel copula is therefore given by: 

\ u = 2 - lim 1 - - _j\ X } = 2 - 25 , and A ; = 



and we have 



g(v,u) 



0+ 1 - if^ 1 (x) 

dC (u, v) 
du 

d exp ( — {—Inuf + (—Inv) 6 



exp 



du 

(—Inu) + (—Inv) 1 



((-Inu) 9 + (-Invf^j 
(—Inu) 



(27) 



Note that (27) is a strictly increasing with respect to v. Its inverse g 1 (a,u) is thus well 
defined. However its inverse <? _1 (a,u) cannot be expressed in an explicit form. Hence we 

cannot derive CoVaRa L ~ l analytically, but we can use in this case we can use numerical 
methods. 

5.3.2. Clayton Copula 

The generator of the bivariate Clayton copula is given by: 

(p(t) = tr e - 1. 

According to proposition 6 we have the following expression for the CoVaRa 



CoVaR s J LZ=l = F- 1 



(o-ifr-l) {Fi(l)y e + 1 
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6. Conclusion 

Managing and regulating the systemic risk is a fundamental problem for financial regu- 
lators and risk managers especially in the context of the current crisis. The must important 
challenge here is the modeling and the quantification of the potential contribution of one 
given individual financial institution to the financial system. One of the main approaches to 
solve this problem is the covar method proposed by Adrian and Brunnermeier in [2] . Where 
the financial system is defined as a portfolio of two items such that the loss of the system 
is represented by a random vector (L 1 L s ) where L % is the loss of the focused financial in- 
stitution i and L s the loss of the financial system s, and the marginal risk contribution of 
the bank i to systemic risk s is quantify by the risk measure ACoVaR s][l which is defined as 

the difference between CoVaRu L ~ VaR ' and CoVaRa ^ ^ (see def. [4]). The problem 

slL^—l 

of the computation and the analysis of the term CoVaR a for a given I is thus very 
important for the implementation of CoVaR especially in the non-Gaussian world, but still 
we do not get any definite solution. As an answer to this problem, we have developed our 
method, based on copula theory, an analytical framework for the implementation of the 

s\L i —l 

CoVaR methodology where the risk measure CoVaRa is expressed in a closed form in 
terms of the marginal distributions F s and F$ separately and the copula C between the fo- 
cused financial institution i and the financial system s. This framework provides an effective 
computation and analysis tool for the systemic risk using CoVaR for a widely used class 
of distribution function and comovement dynamic(see Theorem 13, which captures not only 
linear correlation but also nonlinear tail dependencies between the banks in one financial 
system (which summarise the main features of loss distribution) as opposed to the "linear 
quantile regression" and the formula in |35j where the dependence is modeled only by the 
linear correlation coefficient. In fact our approach allows to analyse the marginal effect of 
Fi, F s and C of the systemic risk. We show for example the systemic risk contribution of i is 
independent of Fi (see ( |TT| )) and highlight in remark [6] some properties of CoVaR according 
to the nature of F s . Our approach can also be used to develop closed formulas for the com- 
putation of related macro-risk measures like CoVaR^id. [2]), ACoVaR^, ACollVaR s \ l 
(cf. [33]), and ACondVaR^ (cf. [33J). 
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